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We study theoretically proximity-induced superconductivity and ferromagnetism on the surface of a topologi- 
cal insulator. In particular, we investigate how the Andreev-bound states are influenced by the interplay between 
these phenomena, taking also into account the possibility of unconventional pairing. We find a qualitative dif- 
ference in the excitation spectrum when comparing spin-singlet and spin-triplet pairing, leading to non-gapped 
excitations in the latter case. The formation of surface-states and their dependence on the magnetization ori- 
entation is investigated, and it is found that these states are Majorana fermions in the d vy -wave case in stark 
contrast to the topologically trivial high-r c cuprates. The signature of such states in the conductance spectra is 
studied, and we also compute the supercurrent which flows on the surface of the topological insulator when a 
Josephson junction is deposited on top of it. It is found that the current exhibits an anomalous current-phase re- 
lation when the region separating the superconducting banks is ferromagnetic, and we also show that in contrast 
to the metallic case the exchange field in such a scenario does not induce 0-71 oscillations in the critical current. 
Similarly to the high-T c cuprates, the presence of zero-energy surface states on the topological surface leads to 
a strong low-temperature enhancement of the critical current. 



PACS numbers: 74.45,+c 

I. INTRODUCTION 

Topological insulatori^E] represent a state of matter which 
is distinct from all other condensed matter systems. Their 
hallmark is the formation of topologically protected, conduct- 
ing edge-states (in 2D) and surface-states (in 3D), whereas 
the bulk retains an insulating behavior (see e.g. Ref. [5] for a 
nice introduction). These surface-states are characterized by a 
topological Z2 symmetry and are robust against disorder and 
perturbations that respect time-reversal symmetry. The key to 
this robustness is the fact that the Brillouin zone of topological 
insulators feature an odd number of Dirac cones (in contrast 
to the even number of such cones in graphene), which ensures 
that backscattering paths due to e.g. impurities always inter- 
fere destructively. 

After their prediction 3 , topological insulators have been ex- 
perimentally observed in HgTe/CdTe quantum wellsP and in 
Bi2Se3/Bi2Te3 crystals 6 . Recently, several works have inves- 
tigated proximity-induced superconducting and ferromagnetic 
order on the surface of a topological insulator 7 12 . It has been 
found that the synthesis of the abovementioned orders in the 
environment offered by the topological insulator yield a num- 
ber of interesting possibilities. On the one hand, it has been 
shown that such hybrid structures can host so-called Majo- 
rana fermions 13 . This class of excitations, in contrast to their 
Dirac equivalent, are their own antiparticles and satisfy non- 
Abelian statistics^ 4 -!. The latter aspect, in similarity to the frac- 
tional quantum Hall effect, has prompted suggestions of us- 
ing such excitations in topological quantum computation ow- 
ing to their robustness towards decoherence effects. It should 
be noted that it recently has been proposed that Majorana ex- 
citations may be generated in semiconductor! superconductor 



heterostructures in the presence of a magnetic field 15 ' 16 ! On 
the other hand, the study of topological insulators also attracts 
interest due to the possibility of unveiling novel transport phe- 
nomena with respect to spin- and charge-transporters 

In RefP, it was recently shown how an interplay between 
unconventional superconductivity and ferromagnetism on the 
surface of a topological insulator would give rise to a number 
of effects with no counterpart in conventional metallic sys- 
tems. For instance, the proximity-effect from a spin-triplet 
superconductor would give rise to gapless excitations in the 
topological insulator since the gap simply renormalized the 
chemical potential. Moreover, it was demonstrated how the 
zero-energy state d 20 ' 21 ' formed due to a cL,-wave order param- 
eter were Majorana fermions in contrast to e.g. the topologi- 
cally trivial high-r c - cuprates, and that the dispersion of these 
states would be highly sensitive to the orientation of a mag- 
netic field. All of these findings demonstrate that qualitatively 
new effects may be expected when superconducting and fer- 
romagnetic order conspire in the environment of a topological 
insulator. 

Motivated by this, in this work we present a comprehen- 
sive treatment of hybrid superconductor! ferromagnet struc- 
tures deposited on top of a topological insulator. In particular, 
we focus on their transport properties and how the Andreev re- 
flection process is altered compared to in conventional metal- 
lic systems. We study both point-contact spectroscopy and 
Josephson junction geometries which are directly experimen- 
tally relevant, and allow for the possibility of unconventional 
superconducting pairing such as p-wave or af-wave. Our pro- 
posed model is shown in Fig. [T] A voltage or current bias 
may be applied to a topological insulator where the surface is 
coated with a ferromagnetic insulator and a superconductor. 
In this way, one may access both conductance spectra and the 
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supercurrent characteristics in order to probe how these are 
influenced by the environment of the topological insulator. 
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The present authors have already published a Letter-^ 1 re- 
porting some of the results in the present paper. In addition 
to the more detailed and comprehensive explanation of the 
derivation and results, there are new results presented in Sec. 



Ill A 4 III B and in the Appendix of this work. 



II. THEORY 

We will employ a Bogolioubov-de Gennes approach to ob- 
tain the bound-states and transport properties of the system 
under consideration. Using a Nambu basis 

^=(¥t.¥;,¥t.¥j) ! (1) 

the Hamiltonian for the surface of a topological insulator un- 
der the influence of either a superconducting or magnetic 
proximity effect reads: 



H 



H (k)+M 



A(k) 



(2) 



-A*(-k) -#o*(-k) -M*J ' 

where we have defined 

flo(k) = v F (ajc x + Oyky) -/j (3) 

and denotes a 2 x 2 matrix. The gap matrix A(k) depends 
on both the orbital- and spin-symmetry of the Cooper pair, 
whereas the ferromagnetic contribution reads M = m • a with 
an exchange field m = (m x ,m y ,m z ). In the superconducting 
region, we set M = while in the ferromagnetic region we set 
A(k) = 0. In the following, we set Vp = 1. It should be noted 
that it is also possible to use a Rashba-like term in the normal- 
state Hamiltonian Ho, i.e. (a x k y — G y k x ) instead of (a x k x + 
Gyky). This gives rise to a different spin-momentum locking 
on the Fermi surface, which was observed experimentally in 
Ref. |22j Nevertheless, all our conclusions below pertaining 
to the excitation spectrum and the qualitative behavior of the 
spin-singlet versus spin-triplet superconducting scenarios are 
independent of whether one uses the Rashba- or Dirac-like 
model for Hq. 

It is instructive to consider how the band-structure on the 
surface of the topological insulator depends on the doping 
level and also the presence of a magnetic field. In the absence 
of proximity-induced ferromagnetism, the energy dispersion 
reads: 



e = ±|k|-( 



(4) 



which is equivalent to the massless, relativistic Dirac fermions 
in graphene. However, we underline once again that the Bril- 
louin zone of graphene contains an even number of such Dirac 
cones whereas that number is odd in a topological insulator. 
We also note in passing that depending on the magnitude of 
the chemical potential fu, it is possible to generate so-called 
specular Andreev reflection 23 on the surface of a topological 
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FIG. 1: (Color online) We consider a topological insulator where su- 
perconductivity and/or magnetic correlations are induced on the sur- 
face via the proximity effect to host materials with the desired prop- 
erties. The possibility of unconventional pairing such as rf-wave is 
taken into account, and the degeneracy of the band-structure is lifted 
in the presence of an exchange field induced from the ferromagnetic 
insulator. 




FIG. 2: (Color online) The band-structure on the surface of the topo- 
logical insulator. The helical surface states reside on a Dirac cone, 
which is influenced in a qualitatively different manner for various 
magnetization orientations. 



insulator between a region with and without a superconduct- 
ing gap Ao since the reflected hole may have parallel group 
velocity and momentum when /j < Ao. 

In the presence of an exchange splitting induced by a mag- 
netization, the band dispersion becomes: 



e = ± 



^(k x + m x ) 2 + (k y + m y ) 2 + 



mj—jj, 



(5) 



It is seen that the exchange field enters in the same way as 
a vector potential would in a conventional metal. While the 
transverse components {m x ,m y } of the magnetization shift the 
position of the Fermi surface, the z-component of the field 
has a qualitatively different effect - it induces a gap in the 
spectrum. This is seen by considering the T point k = 0, where 
the two bands become separated by an energy gap 2m z . The 
influence of the magnetization on the band-structure is shown 
schematically in Fig. [2] As we shall see later, the direction of 
the magnetization strongly influences the transport properties 
of our system. 

We close this section by briefly outlining the relation be- 
tween the Hamiltonian in Eq. (|2|i and the Hamiltonian of the 
bulk topological insulator. The Hamiltonian for a 3D topo- 
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logical insulator such as Bi2Se3/Bi2Te3 can be written down 
by symmetry arguments as shown in RefP2 To account for 
the surface-states, the eigenstates SP are obtained by diagonal- 
izing the full Hamiltonian and then demanding that VP = at 
the boundaries. In this way, one obtains a set of surface-states 
with a linear energy-momentum dispersion which cross at the 
T-point. These surface-states can be influenced by finite-size 
effects 25 27 if the topological insulator has a sufficiently small 
width, typically a few nm in Bi2Se3. In our work, we assume 
that this width is sufficiently large to rule out any finite-size 
effects such that the surface-states simply can be described by 
a Dirac cone, which leads to Eq. Q. 

III. RESULTS & DISCUSSION 
A. Conductance spectroscopy 

In this section, we will focus on the transport properties of 
hybrid structures deposited on top of a topological insulator 
in order to investigate how the interplay between ferromag- 
netism and superconductivity in this context can be probed 
by conductance spectroscopy. To this end, we will employ 
a scattering-matrix formulation along the lines of Blonder- 
Tinkham-Klapwijk 28 theory, including also the role of uncon- 
ventional superconducting pairing. To accommodate super- 
conductivity by means of the proximity-effect experimentally, 
it is necessary to realize the condition /j Aq to have a suf- 
ficiently large density of states. Throughout this paper, we 
will therefore be concerned with precisely this situation and 
set /j S> Ao everywhere on the surface of the topological in- 
sulator except in the ferromagnetic region where we consider 
/j = 0. The reason for this is twofold. Firstly, such a large 
Fermi-vector mismatch effectively accounts for a barrier be- 
tween the non-superconducting and superconducting region 
of the surface of the topological insulator, which is expected 
to be present experimentally. Secondly, it is necessary to use 
an insulating material between the voltage source and the su- 
perconducting region in order to ensure that transport occurs 
exclusively along the surface of the topological insulator (and 
not through the bulk of the host material inducing ferromag- 
netism). In this way, a ferromagnetic insulator deposited on 
top of the topological insulator would be suitable experimen- 
tally, as depicted in Fig. [T] We now proceed to discuss the 
wavefunctions in each of the regions of the surface of the topo- 
logical insulator. 

In the normal (N) region of the surface, we have M = A = 
0, and an incoming right-moving electron with energy 8 may 
then suffer two possible fates upon scattering: i) reflection 
as an electron or ii) Andreev-reflection as a hole. The total 
wavefunction may then be written as: 

y N = e^([l,e* 0,0]e^* + r4l,-e-* 0,0]e^ 

+ r /! [0,0,l,-e- i9 ]e i ^), (6) 

where 8 denotes the angle of incidence while r e and r/, are 
the normal and Andreev scattering coefficients, respectively. 



Here, k x = jUatcosO is the jc-component of the momentum 
which is non-conserved due to the broken translational sym- 
metry, whereas k y = /j^sinG is conserved. 

In the ferromagnetic insulator (F) region of the surface, we 
have A = and M^O. We set the chemical potential to zero 
in this region for reasons outlined above. In that case, the 
wavefunction in the F region is in general a superposition of 
right- and left-moving electrons, due the transmission from 
the N region and reflection at the S interface, in addition to 
right- and left-moving holes, due to Andreev-reflection at the 
S interface and normal reflection at the N interface. Thus, we 
obtain: 

y F = tfry ( ai [ia+ ,1,0, Q] e -^+ im *> 

+a 2 [-ia; 1 ,l,0,0]e^^> 

+a 3 [0,0,ia_,l]e( K -+ im ^ 

+a 4 [0,0, -ioC 1 , l]e- (K - - imr)l ) , (7) 

where aj are the scattering coefficients and we have intro- 
duced: 

k± = yjm 2 + (k y ±m y ) 2 , 

oc± = - [k± - (ky ± m y )]/m z . (8) 

It is instructive to consider in some more detail why 
the wavefunction \\tp has the form of Eq. (FT). Di- 
agonalizing the Hamiltonian in the F region provides 
two eigenvalues for the electron-like quasiparticles £ e ± = 
±yj(k x + m x ) 2 + (k y +m y ) 2 + m 2 , and two eigenvalues e± = 
± J (k x — m x ) 2 + (k y ~ itiy ) 2 + 2 for the hole-like quasiparti- 
cles. Since we are considering an incident electron with e > 
from the N side, where the chemical potential is assumed to 
satisfy /an ^> {Ao,e}, this excitation enters the F region in the 
band ei = \/ (k x + m x ) 2 + (k y + m y ) 2 + m 2 . We illustrate this 
in Fig. pi The Andreev-reflection process occurring at the F|S 
interface gives rise to a hole excitation by removing an elec- 
tron from the band el . The hole dispersion is opposite in sign 
to the electron dispersion of the band it was generated in, and 
therefore in this case follows from e+, which means that its 
group velocity is parallel to its momentum. In effect, this is 
therefore a specular Andreev reflection. By constructing the 
eigenvectors for the excitations in the bands ef + and e?L one 
finally arrives at Eq. |7]l when taking into account both left- 
and right-moving exictations. 

Finally, we write down the wavefunction in the supercon- 
ducting (S) region of the surface which includes a contribution 
from both an electron-like and hole-like quasiparticle. As we 
shall see later, the superconducting spin-triplet case is qualita- 
tively different from the spin-singlet case. The wavefunction 
below is therefore only valid for the spin-singlet case, such as 
s-wave and cZ-wave, and reads: 

Vs = e k y? (t e [e ip , e i (P +e ') , - e i (8'-Y+) , e^e^ 

+t h [l, -e- ie ' , e 1 ^- 9 '-^ , ei{fHr-)] e -^ . (9) 
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FIG. 3: (Color online) Illustration of the scattering processes and the 
bands partaking in these at the N|F interface. An incoming electron 
(1) with e > from the N side can be transmitted (2) to the F region 
in the band e5_. When this electron hits the S interface (not shown), 
it has a finite probability of being Andreev-refiected. In this process, 
an electron is taken (3) from the band z e _ and transmitted to the su- 
perconductor, leaving a hole (4) behind. This hole excitation has a 
band dispersion which is displaced in momentum space with a vec- 
tor (2ra x ,2ra v ) compared to the electron-band. Only the x-direction 
is shown above for clarity and we have set m z = 0. For non-zero m z , 
an energy gap 2m, opens between the Dirac-cones in the F region. 



Note that in writing down the above wavefunction we have 
taken into account the possibility of anisotropic pairing, and 
consequently defined 



m + /m_, u± 



(l± x /e 2 -|A(6')l 2 /e 



A(e ± )/|A(0±)|, e + = e',e_=7r-e'. (io) 



A difference in doping level between the N and S regions is 
accounted for by /j^sinG = /j^sinG', since in an experimen- 
tal situation the S region is often heavily doped /j$ 3> (e,Ao). 
The doping level can be controlled electrically by means of an 
applied gate voltage. If the N region is also doped away from 
the Dirac point, /jn = L>s leads to 6' = 8. In the next section, 
we set fJiy = fJs^ (e,Ao) unless specified otherwise. 

In order to calculate the Andreev-bound state energies and 
the conductance spectra of this junction, we need to solve for 
the scattering coefficients. By matching the wavefunctions at 
each interface, i.e. — \\if at x = and = at x = L, we 
obtain the system of equations Ay = b where we have defined: 



A 



A x J? 2 

y = [r e ,rh,a l ,a2,aT, 1 a A ,t e ,t h ], 
b= [-l } -e ie , 0,0, 0,0,0,0]. 



(11) 



The analytical expressions for the matrices Aj can be found 
in the Appendix. After solving the above system for the un- 
known coefficients y, the normalized conductance G/Gq can 
be calculated according to the formula: 



rn/2 

G= / d9cos9[l 

i-Tt/2 



|r ft (9)| 2 -|r e (e)| 2 ], 



(12) 



and the normalization constant is chosen as Go = G(|eV| ^> 
Aq) as usually done in experiments. We will also be interested 



in the bound-state energies of the junction, which correspond 
to resonant states that persist even in the limit of a vanish- 
ing normal-state conductance. This is modelled by letting the 
width L of the junction become very large compared to all 
other length scales. These resonant states are found analyti- 
cally by identifying the energies e where the probability for 
normal reflection vanishes, i.e. r e = 0. 



1. s-wave singlet pairing 



For a spin-singlet symmetry one finds that 
A(k) = A(k)io v . 



(13) 



The i-wave case has an isotropic order parameter A(k) = Ao, 
and diagonalization of Eq. Q then yields the standard eigen- 
values 



e = W(v f |k| - M 2 + |Ao| 2 , ti = ±1, p = ±1. 



(14) 



Employing the strategy described in the previous section, we 
may calculate the resonant states by looking for energies that 
produce r e = when L^oo. When m y = 0, one is able to 
write down a manageable analytical expression for the proper 
condition related to the formation of bound-states: 



e 2l h++x_ =0. 
Here, we have defined the auxiliary quantities: 
x± = sin 8 + sin(29 + 8) ± [sin 9 + sin(28 + 9)] , 



(15) 



v = (k y -K)/m z , 8 = -iln(v/i), k= \Jk 2 +i 
In the range |e| < Aq, one may write 



p = atan[A (n /l-(e/Ao) 2 /e] ! 



(16) 



(17) 



which upon insertion into Eq. ( fT3j ) yields the following solu- 
tion for the bound-state energy: 



e = A sgn{<:_}/ V / l + C 



where we have introduced: 



C± = tan 



In 



(18) 



(19) 



We have checked analytically that Eq. ( [18) is identical to 
Eq. (7) in Ref. [9] To explore how the magnetization influ- 
ences the bound-state dispersion, we plot in Fig. |4]the bound- 
state energy versus the angle of incidence (or equivalently the 
transverse momentum index) for several choices of m z . It is 
seen that the bound-states have a dispersion only near 9 = 
when the Zeeman field is small, \m z \ <C /J. A zero-energy so- 
lution is seen to be allowed for normal incidence as long as 
m z is finite. It should also be noted that the chirality of the 
bound-states are determined by sgn{m-}. To see this, note 
that to, —> (— to z ) leads to 8 — > 8 + Jt. Due to the symmetry 
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FIG. 4: (Color online) Plot of the bound-state dispersion in the s- 
wave case for several values of \m z \/fi. The chirality of the bound- 
state is controlled by the sign of m z . The solid lines correspond to 
m z > while the dashed lines correspond to m z < 0. a: \m z \/p = 0.1, 
b: \m z \/n = 0.5, c: \m z \jp = 0.9. The superconducting gap Aq is 
used as the fundamental energy unit here, and we set /j/Aq = 100. 



relations x±(8 + 7i) = -%(8) and ln(x_/x+) = -ln(x + /x_} 
follows from the definition of C± that 



s g n {Ct} « sgn{m z }. 



(20) 



To investigate how the presence of such bound-states is 
manifested in the experimentally accessible electrical conduc- 
tance, we plot this quantity in Fig. [5] This is done by solving 
numerically for the scattering coefficients, which enables us 
to consider also the case m y ^ 0. The dependence on the mag- 
netization orientation is shown by considering a pure (a) m x - 
component, (b) m y -component, and (c) m 2 -component. Inter- 
estingly, the dependence in all of these three cases are qualita- 
tively very different, even though the s-wave order parameter 
is rotationally invariant. As we shall see in what follows, this 
stems from the unique band-structure of the surface-states of 
the topological insulator. 

The only common feature all three plots in Fig. [5] have is 
the two coherence peaks at £ = Ao which exist when mj —> 0, 
j G {x 7 y,z}- In Fig. [5ja), the conductance is invariant with 
respect to m x , whereas in (b) it displays a considerable depen- 
dence on the magnitude of m y . To explain this feature, the 
key point to observe is that the fc y -component of the momen- 
tum is conserved in the geometry under consideration since 
translational invariance holds in this direction. In this way, 
the component of the wavevector in the jc-direction of propa- 
gation will in general contain a real and imaginary part, where 
the latter causes the wave to be evanescent (decaying). Even at 
normal incidence = 0, the wave becomes evanescent when- 
ever m y ^ 0. In contrast, m x only influences the real part of 
the wavevector and does not influence the resistance of the 
junction since it is of no consequence for length scale of the 
decaying modes. Therefore, the conductance remains invari- 
ant under a change of m x whereas it is strongly dependent on 
the value of m y . In fact, as m y increases the overall tendency 
of the conductance is that it is suppressed. In (c), we consider 



how the conductance changes when varying the "mass gap" 
component m z . As m z grows, it is seen that the appearance of 
the bound state is manifested by a large enhancement of the 
zero-bias conductance. Therefore, the two finite-energy peaks 
are merged into one zero-energy resonance. Finally, we note 
that the £/ v 2_ v 2-wave case is qualitatively similar to the s-wave 
case. 



2. p-wave triplet pairing 

Turning to the spin-triplet case, the gap matrix now reads: 
A(k) = (d k -o)io,. (21) 

Let us consider a general triplet state dk = A(k)z. Diagonal- 
izing Eq. (ph produces the following eigenvalues: 



e = Tiv f |k|-p x /A< 2 + |A(k)| 



(22) 



This equation is qualitatively different than Eq. ( 14 1 for the 
s-wave case. Namely, the superconducting order parameter 
now renormalizes the chemical potential and the excitations 
remain gapless. From the dispersion Eq. ( [22} follows sev- 
eral anomalous properties. By evaluating the corresponding 
wavefunction, one may conclude that Andreev reflection is 
strongly suppressed at the interface to a non-superconducting 
region since there is no gap in the charge excitation spectrum 
that can retroreflect a hole quasiparticle. Moreover, we have 
verified that for any triplet symmetry the anomalous disper- 
sion Eq. (22 1 is obtained. It also holds even if the Dirac-like 



Hq is replaced with a Rashba-like Hq as mentioned previously. 
Thus, the results for singlet and triplet pairing differ qualita- 
tively in a fundamental way, as the excitations are gapped in 
the former case whereas they remain ungapped in the latter 
case. The structure of the eigenvalue in Eq. ( |22] > appears to be 
a direct result of the band-structure in the topological insula- 
tor, where the spin couples directly to momentum through the 
term o • k in the Hamiltonian. Due to the fact that spin will 
be parallel to the momentum, it follows that pairing between 
equal spins (triplet pairing) at k and — k is not possible. This 
should be distinguished from the case of graphene, where the 
operator a does not represent physic al spin , but rather a pseu- 
dospin index related to the sublattices^ft^U. 



3. d-wave singlet pairing 



We finally address the cZ-wave pairing state, focusing here 
on a d xy -wave symmetry. The reason for why this sym- 
metry is the most interesting is that it is known to produce 
zero-energy surface states in the high-!T c cuprates. The or- 
der parameter may in this case be written as A(k) = A(8) — 
Aocos(20 — Jt/2), and diagonalization of Eq. (j2j) then yields 
the standard eigenvalues 



e = T 1 V(v J? |k|-PA') 2 + |A(e)| 2 ,T 1 = ±l, p = ±i. (23) 
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Proceeding in the same fashion as the previously considered s- 
wave case, we find the following condition for the bound-state 
energies when m y = 0: 



3 2i Px + -T_ =0. 



(24) 



As seen, the only difference from Eq. ( 15 1 is the sign of the 
last term in the equation, although we shall see that this sign 
change has fundamental consequences. Eq. (24i yields the 
following solution for the bound-state energy: 



e=|A(9)|sgn{C + }/Vl + C- 



(25) 



We now consider how the magnetization influences the 
bound-state dispersion in Fig. [6] The dispersion is weak, and 
the bound-state energy is close to £ = 0, for \m z \ <C /J- Increas- 
ing m z (m z > 0) in the of-wave case has the important effect of 
accomodating finite-energy bound-states when moving away 
from normal incidence. It is noted that in the same way as for 
the i-wave case, the chirality of the bound-state is determined 
by the sign of m z . 

The zero-energy states appear in the d xy -wa\e case even in 
the absence of magnetization, as may be shown by solving the 
scattering problem for a N|J AV -wave junction on a topological 
insulator when there is a Fermi-surface mismatch between the 
N and S region. We now show that these zero-energy bound- 
states are Majorana fermions, in contrast to the zero-energy 
states realized in the topologically trivial high-T^ cuprates. 
The essential point in this context is the spin-degeneracy of 
the Fermi surface in the latter case, whereas for a topologi- 
cal insulator this degeneracy is lifted. In both cases, the 4x4 
BdG Hamiltonian H satisfies a particle-hole symmetry 



0//(k)0 = -//*(-k), 
where we have introduced the matrixPH 



© = 



I 

1 



From this property, one may prove that if 

\|/ £ = [wi(k),M2(k),vi(k),v 2 (k)] 



(26) 



(27) 



(28) 



< o 




FIG. 6: (Color online) Plot of the bound-state dispersion in the d xy - 
wave case for several values of \m z \/jj. Similarly to the .s-wave case, 
the chirality of the bound-state is controlled by the sign of m z . The 
solid lines correspond to m z > while the dashed lines correspond to 
m, < 0. a: \m z \/fi = 0.1, b: \m z \//j = 0.5, c: \m z \/n = 0.9. We have 
set ft/ 'Ao = 100. 



is an eigenfunction for the eigenvalue £, then 



© Ve (-k)*= ¥ _ E (k) 

= [vI(-k),vI(-k),«K- 



■k),u2(-k)] (29) 



is an eigenfunction for (—6). For a zero-energy bound state 
£ = 0, one must have \|/ E = \|/_ E , leading to internal symme- 
try relations between the coherence factors such as Ki(k) = 
Vj(— k). The Bogoliubov quasiparticle creation operator for 
this state is constructed in the usual way as 



Y t (k) = Ml (k) Ct t(k)+ M2 (k)c{(k) 
+v 1 (k)c t (-k) + v 2 (k) C; (-k) 



(30) 



Thus, we see that the Majorana criterion y(k) = y^(— k) is sat- 
isfied. Now, the distinction between the zero-energy state in 
the cuprates and the present context of a topological insulator 
is precisely the spin-degeneracy which allows one to split up 
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the 4x4 BdG equations to two separate 2x2 equations per 
spin. Due to the band-structure on the surface of a topologi- 
cal insulator, the e = solution is not spin-degenerate and we 
obtain only one zero-energy mode. As pointed out in RefP, 
this guarantees the Majorana nature of the fermion. We reem- 
phasize that this is different from topologically trivial N| (fry- 
wave junctions, formed e.g. by a normal metal contacted to 
yttrium barium copper oxide (YBCO), where the zero-energy 
solutions are spin-degenerate. 

We now investigate how the presence of such bound- 
states are manifested in an experimentally accessible quan- 
tity, namely the electric conductance. To do so, we consider 
an N|F|<i-wave junction to check how the magnetization can 
be used to manipulate the transport properties and give sig- 
natures of the surface-states. In Fig. [7] we plot the conduc- 
tance, which normally is expected to produce the well-known 
zero-bias conductance peak (ZBCP) due to midgap resonant 
statePEU. 

The common feature for all magnetization directions in (a)- 
(c) is that a zero-bias peak is present when m t ; — > 0, in agree- 
ment with our previous analytical finding. Due to the cou- 
pling between spin and momentum in the band structure of 
the surface of a topological insulator, it is interesting to check 
whether the direction of the magnetization influences the con- 
ductance spectra. In a topologically trivial N|F|4ty-wave junc- 
tion, one can prove analytically that the conductance is invari- 
ant with respect to the direction of the magnetization m of 
the F layer. Increasing the exchange field in the F region, 
the ZBCP splits in the conventional case 33 , similarly to Fig. 
|7|c). We here show that in complete contrast to the topo- 
logically trivial case, the conductance now features a strong 
dependence on the magnetization orientation. We consider a 
magnetization in the x- and y-direction in Fig. |7ja) and (b), 
respectively. It is seen that depending on the magnetization 
orientation, the conductance features three qualitatively dif- 
ferent types of behavior. For m || x, G/Gq is invariant upon 
increasing m x . For m || y, the ZBCP vanishes upon increas- 
ing m y . For m || z, the ZBCP is split upon increasing m z . In 
fact, the evolution of the conductance spectra in (c) is oppo- 
site to the i-wave case upon increasing the magnetization: the 
zero-bias peak is split into two finite-energy resonances. In 
effect, this means that the characteristic features in the con- 



ductance spectra of s-wave and li-wave superconductors can 
be completely reversed by introducing a Zeeman field in the 
topological insulator. The strong sensitivity to the direction 
of m is a new feature compared the topologically trivial case 
which pertains directly to the anomalous band-structure of the 
topological insulator. The difference between the m x and m y 
cases shown in (a) and (b) is explained in the same way as 
for the i-wave pairing scenario. As seen, an increase in m y 
eventually suppresses the influence of superconductivity and 
the conductance is reduced. 



4. Singlet-Triplet mixing 

Due to the lack of inversion symmetry and concomitant 
presence of asymmetric spin-orbit coupling in the topologi- 
cal insulator, one might envision a mixed s + p-wave super- 
conducting state induced by the proximity effect. For a chiral 
p-wave state A(k) = A p e' e (dk-vector along the z-axis), the 
resulting eigenvalues read 

e = r) X /v F |k| 2 + iU 2 + A2 + A2-2pv^, 

R = v 2 F \k\ 2 (iJ 2 +A 2 p )+A^A 2 p cos 2 2Q, r\ =±l, |3 = ±1. 

(31) 

In the limits A, — > and A p — > 0, this reduces to the previous 
expressions in this work. As seen, the presence of an s-wave 
component ensures that the spectrum is gapped, while a cou- 
pling between the s-wave and p-wave components render the 
excitations sensitive to the angle of incidence 0. The ana- 
lytical expression for the corresponding wavefunction in the 
s + /?-wave case is unwieldy, and we defer from any further 
treatment of the transport properties of such a state here. The 
key point we wish to illustrate with this discussion is that the 
s-wave component (in general, the spin-singlet component) is 
necessary to obtain a gapped spectrum, but the interplay be- 
tween spin-singlet and spin-triplet pairing nevertheless gives 
rise to a quasiparticle spectrum which is sensitive to the direc- 
tion of propagation. 
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B. Josephson current 

We now turn to a study of the Josephson current in an S |F|S 
structure deposited on top of the topological insulator. Since 
the F region is assumed to be insulating, such as EuO or EuS 
with band gaps of a few eV, the transport is ensured to take 
place through Andreev bound-states formed on the surface of 
the topological insulator. These bound-state energies can be 
obtained by matching the wavefunctions in a similar way as 
in the previous section. The wavefunction in the F region re- 
mains the same, while we relabel \\f$ —> y's an< ^ concomitantly 
{t e ,th} —> Ki f /j} where the superscript V stands for the right 
region x > L. It then remains to specify the wavefunction in 
the left superconducting region, which reads: 



(a) 



J(P- 



-ikyX 



-4[l,e IB ',-e i (P +9 '^+),e i (P-T+)e^ ) . (32) 



To identify the energy resonances for this system, we look for 
an energy eigenvalue e which gives a non-trivial solution for 
the boundary conditions \\f' s = at x — and = at 
x = L. It can be directly verified by setting up these equations 
that this eigenvalue must satisfy det5Vf = 0, where 



94 



£Wi 94 2 
94 3 94a 



(33) 



and the 2 x 2-matrices 94j are given explicitly in the Ap- 
pendix. After some cumbersome but straight-forward algebra, 
we obtain the following expression for the bound-state energy: 



(34) 



where we have defined 



r(A(|)) = 2cos(A(j) - 2m x L) cos 2 sin 2 5 

- a(cos 2 9 + cos 2 8) [1 - cosh(2KL)] . 
£ = 2 cos 2 8 cos 2 8 — cos 2 — cos 2 8 
+ cosh(2KL) (cos 2 8 - cos 2 0) . 



(35) 



The definition of quantities such as K and 8 are found in Eq. 
([16j, and we have a = — 1 in the s-wave case while a = 1 
in the d^-wave case. Above, we have set m y — since an 
analytical expression for the bound-state becomes unwieldy 
otherwise. The normalized Josephson current may now be 
evaluated according to the standard expression: 



I h = / d0cos0tanh(pe/2) — , (36) 
J-%/2 dA(|) 



r n/2 
it/: 

which upon insertion of Eq. ( [34} produces: 

rV2 A o [g(0)] 2 cos 3 0sin 2 S 



I/I = sin(A(|) - 2m x L) / d0 

J-n/2 



eC[tanh(pe/2)] 



-l 



(37) 
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FIG. 8: (Color online) Plot of the bound-state energy in the Joseph- 
son S|F|S junction with an (a) i-wave and (b) rf tv -wave symmetry. 
We have set ^//A = 100, m z /n = 0.5, and L/% = 0.02. Since m x 
simply corresponds to a shift in A<|), we have set m x = here. 



Above, we have defined A(0) = Aog(0), such that g(Q) = 1 for 
s-wave pairing and§(0) = cos(20 — k/2) for d xy - wave pairing. 
Also, A(j) denotes the superconducting phase difference. Be- 
low, we shall consider both the current-phase relation and the 
dependence of the critical current on the temperature T and 
the junction width L. The critical current is measured experi- 
mentally as 



I c /k = max A( j, (///())• 



(38) 



To make contact with realistic experimental parameters, we 
estimate the Fermi velocity as 2 ^ vf ~ 5 x 10 5 m/s. Due to 
the lattice mismatch between the host proximity materials and 
the topological insulator, the induced superconducting order 
parameter Ao can be expected to be substantially reduced in 
magnitude on the surface of the topological insulator, and 
may be assumed to satisfy Ao ~ 0.5 meV in the li-wave case. 
The superconducting coherence length is then t, ~ 650 nm. 
For a standard s-wave superconductor such as Al or Nb, the 
gap is typically much smaller than in a high-T 6 cuprate, and 
here one might estimate Ao ~ 0.1 meV. To ensure ballistic 
transport, shorter junction are preferable since L must then be 
smaller than the mean free path Z m f p . Also, the formula for the 
Josephson current above is valid under the condition that the 
junction length satifies L cj. We underline that by consid- 



ering Eq. ( 37 i, we may immediately infer that the current no 
longer necessarily vanishes at A(|) = {0,7t} as long as m x ^ 0. 
Therefore, it is possible to actively tune the current-phase rela- 
tion by means of the transverse magnetization component m x 
which enters in the gauge invariant phase-difference between 
the superconducting order parameters. 

In Fig. [8] we consider first the energy of the Andreev-bound 
state in the S|F|S junction, setting for simplicity m x = since 
a non-zero m x simply would correspond to a constant shift of 
A(j). As mentioned previously, we have set m y — for all re- 
sults pertaining to the Josephson current to obtain analytical 
expressions. For normal incidence, it is seen that the bound- 
state energy goes to zero at A(j) = {— k,k} for the s-wave case, 
whereas it is maximal (e = Ao) when — > ±7t/2. This is in 
agreement with the finding of Ref. 9 In the c/ xv -wave case 
shown in Fig. [HJb), the angular dependence of the bound-state 
energy is modified strongly due to the anisotropy of the gap. 
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FIG. 9: (Color online) Plot of the current-phase relation in the 
Josephson S|F|S junction with an (a) .s-wave and (b) rf vv -wave sym- 
metry. We have set /j/Aq = 100, m z /fi = 0.5, and considered the val- 
ues L/% = 0.01,0.02,0.03,0.04 in the direction of the arrow. Since 
m x simply corresponds to a shift in A(j), we have set m x = here. 



(a) (b) 




FIG. 10: (Color online) Plot of the length-dependence of the critical 
current of the Josephson S|F|S junction with T = and an (a) j-wave 
and (b) c/ cv -wave symmetry. We have set /j/Aq = 100 and considered 
m z /n = 0.15,0.40, 1.0 from top to bottom. 



Experimentally, one may probe the Andreev-bound states in- 
directly through their influence on the current-phase relation, 
which we consider in Fig. (|9j. The phase difference can 
be actively manipulated either by means of current-biasing 
the junction or tuning the magnetization component m x . The 
curves are translated to left or right for non-zero m x , depend- 
ing on its sign 37 . As seen, the current-phase relation is qual- 
itatively similar for the s-wave and d xy -wa\e cases, although 
its magnitude is reduced for the latter. This is related to the 
effective weakening of the gap upon Fermi surface averaging 
of its absolute value compared to the s-wave case. It is clear 
from Fig. [9] that the phase-relation is not purely sinusoidal, 
but contains a contribution from higher harmonics. 

Next, we consider the dependence of the critical current on 
the length L of the junction. This quantity is also routinely 
measured in the context of S|F|S junctions. The most no- 
table feature is that there are no 0-71 oscillations in the cur- 
rent, in spite of the presence of an exchange field m z in the F 
region. To explain this, one should note that the influence of 
the exchange field m z is fundamentally different in the present 
scenario where we consider states residing on the surface of 
a topological insulator as compared to a normal metal. In 
the present case it induces a gap in the spectrum whereas in 
the latter case it splits the energy-bands of the majority- and 
minority-spins. As a result, the exchange field does not in- 
duce any finite center-of-mass momentum for the Cooper-pair 
and hence a monotonous decay of the critical current versus 
L should be observed regardless of the magnitude of m. The 
main difference between the s-wave and c/ ry -wave case in Fig. 



10 is that the current becomes suppressed more rapidly as L 
grows in the latter scenario. 

It is natural to next address the question: what happens 
when the magnitude of the exchange field m z decreases and 
eventually vanishes? This would in effect render the junction 
into anS|N|S system. Several previous works 34 35 have inves- 
tigated a similar scenario when using a <f -wave model relevant 
for the high-r c cuprates. In such a case, it was shown that 
the critical current becomes strongly enhanced at low tem- 



peratures in the J ly -wave case, whereas it saturates in the 
s-wave case. The reason for this is the existence of zero- 
energy Andreev-levels in the d xy -wa\e junction. Above, we 
considered the case m z ^> Ao where such states are shifted 
to finite energies as indicated by our previous results for the 
conductance spectra. Now, we study how the temperature- 
dependence of the critical current evolves when m z decreases 
and compare the s-wave and c( xy -wave scenarios. This is 



shown in Fig. 11 



10 the current now 



In contrast to Fig. 
behaves qualitatively different in (a) and (b) corresponding to 
the s-wave and ii„-wave case. In (a), the current saturates 
at a constant value as T/T c —> 0, whereas in (b) there is a 
strong enhancement in the same regime which is more pro- 
nounced the smaller m z /fj becomes. Thus, such an anomalous 
temperature-dependence serves as a signature for the zero- 
energy states on the topological surface in the same way as it 
does for the cuprates, and can be probed experimentally. We 
have verified that when m z 3> Ao, the temperature-dependence 
of the current is qualitatively the same in the s-wave and d xy - 
wave cases. 

We note that an inclusion of the orbital effect due to the 
vector potential A simply would add a component to the mag- 
netization vector as a result of the linear energy-momentum 
dispersion. The predicted results in this work can be tested 
experimentally by fabricating a hybrid structure such as the 
one shown in Fig. [T] In terms of actual materials, EuO 
or EuS might be suitable as ferromagnetic insulators in this 
context 36 . For the li-wave superconductor, a high-T c cuprate 
such as YBCO would be appropriate. Finally, we also point 
out that the suppression of the superconducting order param- 
eter near the interface region has not been taken into account 
here. Such an approximation is valid when there is a strong 
Fermi-surface mismatch, as considered throughout this paper. 
Nevertheless, it could be interesting to see if the gap suppres- 
sion is able to host additional bound-states near the interface 
by employing a self-consistent solution of the order parame- 
ter. We leave such issues for future work. 
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(a) (b) 




T/T c T/T c 

FIG. 1 1 : (Color online) Plot of the temperature-dependence of the 
critical current of the Josephson S|F|S junction with L/% = 0.03 and 
an (a) s-wave and (b) d vy -wave symmetry. We have set ///An = 100 
and considered m z /jj = 0.02,0.08,0.14 from top to bottom. As 
in z /jj —¥ 0, the junction becomes an effective S|N|S system. 

IV. CONCLUSION 

In summary, we have considered the interplay between 
magnetic order and unconventional superconducting pairing 



on the surface of a topological insulator. We find that the 
charge excitation spectrum is rendered gapless for any spin- 
triplet state, such that bound-states are absent and Andreev 
reflection is strongly suppressed. For spin-singlet pairing, we 
find that the zero-energy surface states in the of TV -wave case are 
now Majorana fermions, in contrast to the case of the topolog- 
ical^ trivial high-!T f cuprates. We have studied how Andreev- 
bound states and Majorana fermions are influenced by the in- 
ternal phase of the superconducting order parameter, and find 
that the ZBCP being the hallmark of the d^-wave state is qual- 
itatively strongly modified in the present context. In particu- 
lar, it is highly sensitive to the magnetization orientation, in 
contrast to the topologically trivial case. Our findings can be 
directly tested through tunneling spectroscopy measurements, 
and we have estimated the magnitude of the necessary exper- 
imental quantities. 
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Appendix A: Matrices J?; and Mj 



Introducing the quantity 8± = — iln(oc±/i), we may write 
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which are used to calculate the bound-states and conductance spectra in the N|F|S junction. To obtain the Josephson current in 
the SlFlS case, we use: 



9vt\ 
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(A2) 



Note that in writing down these matrices, we have explicitly separated out the superconducting U(l) phase corresponding to the 
broken symmetry and defined — ioc± = e l8± . Without loss of generality, we set it to ±<|>/2 on the right and left side, respectively. 

The phase-factors y± thus only contain information about the internal phase in k-space for the order parameter. 
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